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Змістовні математичні постановки конкретних геофізичних задач загалом зводяться до визначення 
коефіцієнтів відповідних рівнянь математичної фізики. Ця ідея використана, щоб дедуктивно звести 
нелінійне інтегральне рівняння у постановці задачі Алексідзе до лінійної комбінації шуканих розв'язків. У 
такій альтернативній постановці задача зведена до класичної задачі варіаційного обчислення. 


Содержательнье математические постановки конкретньх геофизических задач в общем сводятся 
к определению козффициентов соответствующих уравнений математической физики. Зта идея испо- 
льзована, чтобьї дедуктивно свести нелинейное интегральное уравнение в постановке задачи Алексидзе 
к линейной комбинаций искомьх решений. В подобной альтернативной постановке задача сведена к кла- 
ссической задаче вариационного исчисления. 
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Обмеження математичного моделювання в геофізиці. Досі представники функці- 
онально-аналітичного напрямку геофізики розробляють методи і алгоритми у рамках 
відомих інтегро-диференціальних рівнянь математичної фізики, які описують поведін- 
ку геофізичних полів. Їхні зусилля спрямовані на те, щоб визначити деякі коефіцієнти 
цих рівнянь та внести відповідні поправки за "складність" середовища. Ці коефіцієнти 
визначаються за допомогою ітераційних схем, які вкладаються у принципи Тихонова 
розв'язування обернених задач математичної фізики (4). 

Одначе, визначення цих коефіцієнтів за функціями, які утворюють відповідні систе- 
ми рівнянь - принципово неоднозначне. Це особливо очевидно при масових машинних 
розрахунках на системах великої розмірності |7|. А за найменших похибок у вхідних 
даних розв'язки - істотно нестійкі. Подолати цю нестійкість на порядок важче, ніж от- 
римати формальні розв'язки задач. До того ж, слід вхідних функцій, за яким визначають 
їхні коефіцієнти - двовимірний, а відновити необхідно тривимірне поле, за яким знахо- 
дять коефіцієнти при третій координаті, яка не співпадає з розмірністю сліду. 

Формально підійти до цієї проблеми означає отримати псевдопросторовий розподіл, 
який "розчиниться" серед еквівалентних розв'язків задачі. Крім того, запозичивши тео- 
ретичний апарат математичної фізики, геофізики одночасно успадкували суто евклідову 
структуру просторів функцій, яка для регіональних побудов втрачає свої якості |21. 
Отже, обмеженість інтегро-диференціальних операторів, які апроксимують математич- 
ну модель середовища - це методичний зріз проблеми 

Вхідні дані - це, по-суті, двовимірні проекції геофізичних полів на поверхню Землі, 
виміряні з похибками на розрідженій нерегулярній мережі спостережень. Вони апріорі 
отримані з невідомими похибками, тому додатково до сказаних вище методичних проб- 
лем ускладнюють поведінку (1 визначення) функцій, які входять до відповідних модель- 
них рівнянь. А розвиток комп'ютерного моделювання призвів до загального застосува- 
ння дискретної шкали вхідних функцій. Зважаючи на те, що досі не освоєно пряме от- 
римання вхідних даних для обернених задач геофізики з баз даних пунктів спостере- 
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жень |3), переведення неперервних вхідних даних у дискретну форму в процесі цифро- 
вої картографії неминуче вносить додаткові спотворення, окрім похибок вимірювань. 
Ці проблеми є наслідком методології минулого століття, що відживає своє, тому ми їх 
називаємо методологічними. 

Новий геофізико-математичний апарат. Як вказані вище, так і інші математичні 
та методичні проблеми, зумовлені неадекватністю апарату математичного моделюван- 
ня геофізичних даних, змусили геофізиків розробляти власний геофізико-математичний 
апарат. В його рамках істотно перероблені відомі математичні методи з урахуванням 
фактів, притаманних суто геофізичним явищам. Започаткував розробку своєрідного 
"геофізичного діалекту" В.М. Страхов |6| - у потенціальних полях, врахувавши орто- 
гональність сигналу й похибки та дискретність даних; продовжив В.В. Гольдін - у сейс- 
міці, врахувавши нелінійність акустичного середовища та сейсмічну анізотропію кори; 
підтримав В.В Аксьонов |1| - у електродинаміці, врахувавши сферичність джерел елек- 
тромагнітного поля і механізми динамозбудження, 1 т.д. 

Без чіткої національної програми розвитку технологій потенціальних полів та ство- 
рення мережі моніторингу гравімагнітних полів ми не можемо вплинути на реформа- 
цію принципів отримання вхідних даних. Сподіваємося, що адекватна потребам як гео- 
фізиків-польовиків, так і фахівців-інтерпретаторів система неперервних спостережень 
для формування загальнодоступної бази даних гравімагнітних полів постане в найшли- 
жчій перспективі. Але наразі ми можемо запропонувати деякі міркування щодо нових 
підходів до зображення математичних моделей геологічного середовища та геофізик- 
них полів. Наші пропозиції вкладаються в рамки нового апроксимаційного підходу |3| 
до створення вказаних конструкцій. 

В рамках дослідження |3| отримано нову аналітичну (диференціальну) модель поля 
сили тяжіння. Основні етапи її виведення ілюструє такий ланцюжок рівнянь: 
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В рамках розв'язання однорідної частини цього рівняння постала так звана задача 
Алексідзе для рівняння Лапласа. Нагадаємо, що для визначення модуля градієнта сили 
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а для визначення поверхневої густини б с є ОС, яка входить у вираз (2), таке нелі- 
нійне рівняння: 
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де спрямовуючі косинуси компонент градієнта елегантно визначаються через суму нор- 
мованих скалярних добутків векторів нормального і аномального поля, тобто со8 р, 4) - 
це у -є) ба -л,) 

га | б | 719 

Методика розв'язування. Для отримання розв'язку рівняння (3) в |3| окреслено па- 
ру ітераційних схем, але вони мало придатні для практичних обчислень. На нашу 
думку, це відбувається в силу слабкої збіжності відповідного інтегралу (в смислі голов- 
ного значення). А його регуляризація через зображення у вигляді відрізка ряду не гара- 
нтує бажаної точності розв'язку. 

Один із можливих способів розв'язання рівняння (3) може бути таким. Зобразимо 
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Тоді його квадрат с с ( )- У ага ех х)ди(х ). Тепер спробуємо формалізувати вираз (3) 
через ее позначення: 
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У такому зображенні задача у підсумку зводиться до класичної варіаційної постанов- 
ки: визначити коефіцієнти а, рівняння (4) таким чином, щоб виконувалась умова міні- 


мізації нев'язки 
2 . 
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Відтепер, залежно від того, яку з популярних у теорії обернених задач метричних 
норм ми оберемо, задача знаходження аналітичної апроксимації (4)-(5) розв'язується у 
тому чи іншому просторі за допомогою класичних варіаційних методів, зокрема, моди- 
фікацій методу Ньютона. 
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